

install.packages("strucchange")

library(strucchange)

#setwd("~/Dropbox/Job Market Paper/IO Revisions/Replication_Files/Replication_Appendix/A_Fig_3")


data_q<-read.csv(file="data_fig_A3.csv")




####### Figure A3 #########


num_states<-tapply(data_q$ID,data_q$Year,length)

years<-as.numeric(as.character(unique(data_q$Year)))




par(mar=c(4,4,1,4))


c2<-Fstats(num_states~1,from=.05,to=.9)



pdf(file="Figure_A3.pdf")

plot(c2,axes=F,xlab="Year",xlim=c(0,1),lwd=2)

axis(2)

axis(1,at=seq(0,1,by=1/7), label=seq(1100,1800,by=100)    )

par(new=T)

plot(years,num_states,type="l",lty=2,axes=F,xlab="",ylab="",ylim=c(100,275))

axis(4,at=seq(100,300,by=25))


dev.off()
